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Probing the evolution of physical systems at 
the femto- or attosecond timescale with light re- 
quires accurate characterization of ultrashort op- 
tical pulses. The time profiles of such pulses 
are usually retrieved by methods utilizing opti- 
cal nonlinearities, which require significant signal 
powers and operate in a limited spectral ranged. 
We present a linear self-referencing characteri- 
zation technique based on time domain localiza- 
tion of the pulse spectral components, operated 
in the single-photon regime. Accurate timing of 
the spectral slices is achieved with standard single 
photon detectors, rendering the technique appli- 
cable in any spectral range from near infrared to 
deep UV. Using detection electronics with about 
70 ps response, we retrieve the temporal profile of 
a picowatt pulse train with ~ 10 fs resolution, set- 
ting a new scale of sensitivity in ultrashort pulse 
characterization. 

Complete characterization of an optical pulse implies 
the knowledge of its time profile, or, equivalently, its 
spectrum and the spectral phase. Despite the rapid 
progress of detection electronics^, direct measurement of 
a subpicosecond time profile remains a difficult task. The 
spectral phase can be reconstructed from the interference 
with a known coherent reference pulsed. Usually, how- 
ever, such reference is unavailable, and the pulse shape 
retrieval requires a self-referencing technique. 

Self-referenced characterization of an ultrashort pulse 
is usually performed by interacting it with its replica in 
a nonlinear optical medium. During the last two decades 
nonlinear optical characterization techniques such as 
FROG 5 and SPIDER 6 have become a standard tool in ul- 
trafast research. Employing optical nonlinearities, how- 
ever, imposes substantial limitations on power, spectral 
range and complexity of the characterized pulsed. 

These difficulties stimulated attempts to develop spec- 
tral phase retrieval methods using only linear optics 
and "slow" detectors. The ability to characterize pulses 
with a temporal resolution significantly better than the 
detector response is intimately related, as in preci- 
sion spectroscopy, to the system noise characteristics. 
Thus, while there are several approaches that allow for 
reliable picosecond waveform reconstruction^ - — , few 
works have demonstrated linear optical characteriza- 
tion with femtosecond resolution, whether by spectral 
shearing^, time resolved^ interferometry, or direct*^ 
or cross correlation^ sonogram measurements. In all 
these, the extracted quantity is the derivative of the 



spectral phase <p(u>) with respect to frequency u. This 
can be thought of as the arrival time of a spectral slice 
of a pulse, provided that the slice is narrow enough so 
that the phase within it can be approximated by a lin- 
ear function^ 5 .. The best reported temporal resolution 
in linear chracterization was of 50fs at a power level of 
O.lmwattii^. All linear methods significantly fall behind 
the few nanowatt sensitivity recently demonstrated in 
nonlinear characterization^. Yet, as we show below, this 
is not a true limitation of linear characterization. 

In this Letter we demonstrate a linear self-referencing 
ultrashort pulse characterization technique with shot 
noise limited temporal resolution. The technique relies 
on precise timing of spectral slices using photon counting 
instrumentation. We take advantage of the fact that, de- 
pending on the signal to noise ratio, the timing error can 
be much smaller than the detector response time. For 
a train of identical pulses, the signal to noise ratio can 
be improved by averaging over many pulses. The overall 
resolution then increases as V^V, where N is the number 
of photons detected (shot noise resolution scaling) . Uti- 
lizing a differential measurement scheme with the signal 
and the reference pulses coupled into the same detector, 
we demonstrate the v^V uncertainty scaling down to few 
femtoseconds. The principle of converting timing accu- 
racy into temporal resolution is analogous to superresolu- 
tion microscopy techniques relying on spatial localization 
of individual fluorophores^— . 

We implemented this method using a fast single pho- 
ton avalanche photodiode (idQuantiquc idlOO) for signal 
detection and a time correlated single photon counting 
(TCSPC) module (PicoHarp 300) for electric pulse tim- 
ing. The ability of TCSPC electronics to detect delays as 
large as few nanoseconds allowed for very large temporal 
dynamic range. The layout of our setup is presented in 
Fig. [U The FWHM instrument response was measured 
to be about 70 ps. Spectral slices were selected by a 4f 
pulse shaper— with a mobile slit in the Fourier plane. 
At every slit position, the data was collected for 0.5s; 
the spectrum was scanned several times to achieve the 
desired integration time. For every slice, we generated 
a histogram of arrival times of the signal relative to the 
trigger pulse. The arrival time was obtained by finding 
the centroid of the histogram peak (see Supplementary 
Section 1 for details of the data treatment). 

The femtosecond relative delays between the spectral 
components are much smaller than the systematic shifts 
induced by the intensity dependence of the detector re- 
sponse and the drifts of the centroid due to temperature 
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FIG. 1: a) Experimental setup. The characterized pulse 
is passed through a 4f pulse shaper selecting a spectral com- 
ponent, and is detected by a single photon counting APD. 
One replica of the original pulse is used as a trigger signal 
for the TCSPC module, while the other, serving as an addi- 
tional reference, is directed into the signal detector. The inset 
shows a sample histogram with the two peaks corresponding 
to the signal and the reference pulses, b) Allan variance. 
The plot shows the standard deviation of the measured pulse 
arrival time as a function of the integration time. Two repli- 
cas of the same pulse separated by a constant delay were 
used as the signal and the reference, each attenuated to 1 
MHz detection rate. The arrival time of the signal pulse was 
measured by finding the absolute position of the peak in the 
histogram (magenta) and by finding its position relative to 
the reference pulse (blue). The dashed black line represents 
the Cramer- Rao minimum variance for this measurement (see 
Supplementary Sections 1 and 2) 



fluctuations. To overcome this difficulty, we introduced 
another reference beam: a delayed replica of the original 
pulse was directed into the detector along with the signal 
beam. Since the photon flux was lower than the satura- 
tion threshold, the detector was randomly activated by 
photons coming from one of the two beams. The ar- 
rival time of each spectral component was then given by 
the relative delay of the two corresponding peaks in the 
histogram. An Allan variance plot, demonstrating the 
measurement timing error as a function of the integra- 
tion time for a single peak and for the interval between 
the two peaks is shown in Fig[T] b. While the unrefer- 
enced measurement deviates from the y/N scaling by an 
order of magnitude already at 1 s integration time, the 
relative delay measurement remains shot noise limited 
up to ~ 1000 seconds of integration, corresponding to a 
sub-cycle standard deviation of 2 fs. 

We first tested our setup by characterizing test pulses 
from a Ti:Sapphire mode-locked oscillator. We retrieved 
the 85 fs transform limited laser output, as well as a 
chirped pulse obtained by passing it through a 6" slab 
of F3 glass. For comparison, the same pulses were char- 
acterized by a second harmonic FROG setup utilizing a 
100 /im BBO crystal. For the linear measurement, the 
characterized beam was attenuated to about 1 photon 
per pulse (approximately 10 pW), leading to ~ 10 6 de- 
tector counts per second. The integration time was 40 s 
per data point. Prior to the measurements, a calibration 
curve accounting for any uncompensated dispersion in 
the setup was recorded using a transform limited pulse. 

The measured spectral component delays with and 
without the glass slab are presented in Fig. [3J The res- 
olution of the measurements, estimated by calculating 
the standard deviation in the range 784 to 789 nm after 
subtraction of the linear trend, was 9.3 fs. for the glass 
measurement and 10.6 fs for the transform limited pulse. 
These values are in good agreement with the y/N scal- 
ing of resolution discussed above (see also Supplementary 
Section 1). As expected, the plots show linear delay de- 
pendence on wavelength, with the slope of — 41.14fs/nm. 
The 10% systematic deviation from the results of the 
FROG measurements (— 46.4fs/nm) and from Sellmeier 
equation for the F3 glass^, (— 47.4fs/nm), which are also 
shown in Fig. [2] can be attributed to a dependence of 
the calibration curve on the exact alignment of the beam 
into the detector. This coupling sensitivity is probably 
due to the fact that in our detector the entrance aper- 
ture is connected to the actual APD by a short length 
of a multimode fibre. The temporal pulse profiles recon- 
structed from our measurement and from the FROG data 
are shown in Supplementary Fig. 2. 

The system's ability to reconstruct complicated pulses 
was tested by passing the beam through a hot Rb vapor 
cell. Rubidium has two narrow absorption lines in the 
relevant spectral range, at 780 and 795 nm. At large op- 
tical densities, achieved by heating the cell, the transmit- 
ted pulses are distorted by the resonant lines dispersion 
profile. The resulting complex pulses contain temporally 



3 



CD 



4 







x 10 s 




10 






8 5 












4100 4200 4300 
time (ps) 


1 







775 780 785 790 795 800 



600 F* s ; * . 




775 780 785 790 795 800 

wavelength (nm) 



FIG. 2: Characterization of a 85 fs transform limited 
pulse and a pulse chirped by a 6" long F3 glass slab, a) 

Density plot of a typical raw measurement trace featuring two 
peaks, the reference and the signal. The inset shows a cross 
section of the signal peak with FWHM of about 70 ps. b) 
Measured delays of the spectral components for the transform 
limited and the chirped pulse (represented by the blue and red 
dots, respectively). The dashed black and cyan lines show the 
delays derived from the FROG measurements, and the dashed 
green line represents Sellmeier equation. The gray solid curve 
in the lower part of the plot shows the spectrum of the pulses 
in arbitrary units. 



extended tails, which make their nonlinear optical char- 
acterization challenging. Each resonance modifies the 
electric field of the pulse multiplying it by a factor of 
f(uj) = exp [— al (1 /(l + 11JJT2) )], where T2 is the inverse 
pressure broadened linewidth, al is the optical density 
and uj is the frequency detuning^. The corresponding 

spectral slice delay is then t(uj) — alT2 ■ The 

results of the measurements are presented in Fig. [3] The 
spectral resolution of the system was set to 0.5 nm, thus 
the observed relative delays were weighted averages of the 
true delays over that range. The spectral delay measure- 
ments were complemented by a higher resolution absorb- 
tion spectrum measurements shown in Fi^SJ:. At room 
temperature, no effect is observed. At 100°C the delay 
peaks around the two resonances are clearly visible, al- 
though the shape of the delay peaks is dominated by 
the instrument response. As the temperature of the cell 
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FIG. 3: Measurement of the spectral delay of pulses 
passed through a Rb cell, a) Measured delays of spectral 
components at 100°C (blue) and 210°C (red). The integration 
time was 4s per point at 100°C and 8s per point at 210°C. 
b) A magnified view of the same (markers) fitted with a res- 
onance delay profile (lines), c) Normalized spectrum of the 
transmitted pulse at 100°C (blue) and 210°C (red) measured 
at a resolution of 0.05 nm. (The dashed black line shows the 
resonant absorbtion profile corresponding to the fit of 210°C 
spectral delay data.) 

is increased to 210°C, the delay in the Lorentzian tails 
of the above profile becomes significant and is clearly ob- 
served. According to the equations above, the magnitude 
of the tails of the spectral delay curve is aL (T2UJ 2 ) , 
while the magnitude of the Lorentzian tails of absorb- 
tion lines in the transmitted spectrum (shown in Fig 
|3fc) can be written as 2aL (T2w)~ . Extracting the two 
coefficients for the 780 nm resonance from the 210°C 
data (see Supplementary Section 4) and comparing them 
gives an estimate of self-broadened resonance linewidth 
ofr = l/X2~10 GHz, which is in reasonable agreement 
with the spectroscopic data on Rb pressure broadening 24 . 
The maximal delays measured in this experiment ex- 
ceeded 30 ps, corresponding to a temporal dynamic range 
of about 1000. 

While the above experiments demonstrate the possi- 
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bility of characterization of ultrashort pulses with 10 fs 
resolution with just a few picowatts of signal power, the 
performance of the method is not limited by these val- 
ues. The maximal resolution and sensitivity achievable 
in this approach are determined by the integration time, 
which, in turn, is limited by the stability of the charac- 
terization setup and of the source itself. This potentially 
enables complete characterization of ultrashort scatter- 
ing from microscopic sources, such as single molecules or 
nanoparticles, a regime relevant for coherent nonlinear 
microspectroscop y 25 i 26 . 

The integration time necessary to achieve a given level 
of temporal resolution can be reduced by using an ar- 
ray of fast detectors, either proportional or single pho- 
ton counting, for spectrally multiplexed measurements. 
Another possible direction that can be explored is the 
characterization of ultrashort pulses in the UV and XUV 



range. The accumulative nature of this method makes 
it potentially suitable for characterization of attosecond 
pulse sources with megahertz repetition rates 27 . 

In conclusion, we presented and demonstrated experi- 
mentally a self-referencing ultrashort pulse characteriza- 
tion technique that combines the single-photon sensitiv- 
ity and versatility of linear detection with the temporal 
resolution previously achieved only by nonlinear optical 
methods. This technique can become a potent tool for ul- 
trafast light sources characterization and for the analysis 
of the fast dynamics of microscopic physical systems. 
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